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Abstract 

We generalize the method of obtaining the fundamental linear partial dif- 
ferential equations such as the diffusion and Schrddinger equation, Dirac and 
telegrapher's equation from a simple stochastic consideration to arrive at cer- 
tain nonlinear form of these equations. The group classification through one 
parameter group of transformation for two of these equations is also carried 
out. 
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I. INTRODUCTION 



It is remarkable that some of the fundamental linear equations of physics such as the 
diffusion and Schrodinger equations, Dirac and telegrapher's equations, and the Maxwell 
equations can be obtained by setting up a master equation from simple stochastic consider- 
ations and a modification there of 

Hosts of nonlinear partial differential equations have been vigorously studied in recent 
times for possessing interesting solitonic, self focussing and allied properties in fields wide 
ranging from lasers to string theories. It is interesting to note, as it is shown in classical 
mechanics with randomness playing a key role, that any arbitrary initial distribution of 
velocities must eventually become Maxwellian. Thus the Maxwellian distribution must be 
invariant under the underlying stochastic process ||. The system tending to and remaining 
in Maxwellian distribution is analogous to the self focussing in nonlinear optics and the 
properties of solitonic systems arising in nonlinear phenomena. 

The transition to nonlinear equation from a linear equation through certain transforma- 
tions is well known. An example is the Cole-Hopf transformation that carries over the linear 
diffusion equation to the nonlinear Burgers equation. The reverse process of getting linear 
diffusion equation from a nonlinear diffusion equation in the form ^ = f 2 §^t has also been 
studied through a nonlinear transformation. 

In this paper, we propose to obtain a class of nonlinear equations in a different way by 
generalizing the method of 0. The method is simple. The basic inputs can be incorporated 
from considerations and arguments based on physical reasoning to obtain nonlinear equations 
rather than arbitrary mathematical transformations. The form of the equations obtained 
are quite restrictive. However, we do not address the deep mathematical significance of 
Cole-Hopf transformation and the like in this method. 

In section II, we briefly review the method used in setting up some basic linear equations 
of physics. Then we generalize the procedure of obtaining classes of correspondng nonlinear 
partial differential equations in section III. Next, section IV is devoted to the construction 
of the groups under which two of the equations obtained in section II, namely diffusion 
equation with nonlinearity and nonlinear telegrapher's equation, remain invariant. The 
similarity transformation and the Lie algebra are constructed to show the transformations 
under which solutions go over to new solutions. 



II. LINEAR EQUATIONS 



The fundamental role of linear diffusion equation in Physics and its applications in other 
branches hardly needs any elaboration. An analytic continuation turns this equation into 
Schrodinger equation. Another basic equation of physics, the Dirac equation originated in 
the attempt to make Schrodinger equation compatible with special theory of relativity. It 
is a simple matter to arrive at telegrapher's equation by iterating Dirac equation in (1 + 1) 
dimensions. It is a curiuos fact that these equations have been obtained by Gaveau etal 
and Ord ||] from stochastic consideration by setting up a master equation. Ord || has also 
arrived at the Maxwell equation by a modification of the master equation. 

Following [1|-|3| we briefly review how these equations are achieved and then proceed to 
nonlinear generalization. The basic consideration is the correlation over a random ensemble 
of particles. However, for simpler visualization we may follow the Boltzmann approach by 
analysing the movement of a single particle. Let a particle have random motion in one space 
dimension moving with a fixed speed v. We assume that it has complete reversal of direction 
of motion in a random manner from time to time, say as to flip of a coin. So this is according 
to Poisson distribution, that is to say that there is a fixed rate a for this reversal and the 
probability for reversal in a time interval dt is adt. Let P+(x, t) (respectively, P-(x, t)) be the 
probability density for the particle being at x at time t and moving to the right (respectively 
left). The master equation for an infinitesimal time step is 

P ± (x,t + At) = P±(x =F Ax,t)(l -aAt) + P T (x±Ax,t)aAt (1) 

This equation gives rise to the linear equations such as Dirac, telegrapher's, diffusion or 
Schrodinger equations in the lowest approximation under various circumstances. 
To the lowest order in Ax and At, equation (1) gives, 

dP± frt „ x dP± Ax 

and the telegrapher's equation follows by iteration, 

d 2 P± 2 d 2 P± dP ± 

The one dimensional Dirac equation is obtained from (1) by analytic continuation. First 
we identify P± with u±, v <-»• c, % -^- a, and then perform a phase transformation u[x, t) = 

intent , * _ . _ 

e s w{x,t). 1 his results m 



z/i— — = mc c^W — icna z — — . (4) 

In the Feynman path integral formulation through checkers moves on space-time, 1 has 
to be replaced by a factor 1 + ( H ^) At for each step on which a reversal does not take place, 

2 

whereas for reversals there is a factor —iAt 1 ^-. 

The Dirac equation in (1 + 1) dimensions, having two components, has the similar time 
and space dependence in this stochastic approach. But for a single component object we 
find that the linear diffusion equation results, which shows the asymmetry in derivatives 
arising out of the random walk problem. 

A generalization to three space dimensions has been carried out in [|IJ. 

McKeon and Ord || have shown that if movements backward and forward in time is 
superposed as well on the previous motion, then Dirac equation in one dimension results 
without recourse to direct analytic continuation. 

To arrive at the linear diffusion equation in a simple way, we put P± = P T = P, and 
a = 2^7. The master equation (1) reduces to 

P(x, t + At) = ^P(x - Ax, t) + ^P(x + Ax, t). (5) 

Expanding this in a Taylor series about the point (x,t) gives, 

„, , dP(x,t) A n/ , d 2 P(x,t) (Ax) 2 

P(x,t) + —^At + ... = P(x,t)+ 2 } +... (6) 

and equating the lowest order terms we get, 

8P_ _ d 2 P (Ax) 2 _ d^P_ 

~dt ~ dx 2 2At a? U 

where, D = 

It may be noted that the above equation in the context of Brownian motion can be 
obtained from a consideration of one dimensional random walk with a Bernoulli distribution 
of probability and the statistical considerations sets Q 

lim ^—t^— = D (8) 
At^o 2At v ' 

where D is a constant. 

A formal analytic continuation (e. g. t — > it, or D — > ih) leads to Schrodinger equation 
for free particles. A potential V(x,t) can be included by adding a term V(x,t)P(x,t)At to 
the right hand side of equation (5). 



Ord P| has obtained the Maxwell equations in 1+1 dimensions by a modification of the 
master equation. We follow his procedure to show how it is done. First equation (1) is 
modified to, 



P±(x, t + At) = P±(x =F Ax, t) + a(x, t)At (9) 

where, a(x, t) is interpreted as a source and linear combinations of P + and P_ will correspond 
the potentials A(x,t) and $(x,t). To the lowest order in Ax and At equation (9) gives, 

dP±(x,t) A dP ± A , , A 

'-At = +— - Ax + a(x, t)At. (10) 



to 



Writing, 



dt dx 



A(x,t) = -[P+{x,t) + P-(x,t)] 



equation (10) implies, 



$(x,t) = ±[P + (x,t)-P-(x,t)] (11) 



dA(x,t) d$(x,t) . s 

-A_a = - c ^l + a(x,t) (12) 

d$(x,t) dA(x,t) 

= - C ^x~ (13) 

where we have put ^ = c. 

Equations (10) and (11) may be decoupled by differentiating the first with respect to t 
and the second with respect to x to give , 

d 2 A(x,t) _ 2 d 2 A(x,t) da(x,t) 

~~dt 2 ~° dx 2 + dt { ' 

and similarly we get, 

a 2 $(x,t) _ 2 9 2 $(x,t) da(x,t) 

dt 2 ~ ~~dx~ 2 c ~dT~ ( } 

Equations (13), (14), and (15) are equivalent to Maxwell equations in (1 + 1) dimensions, 
equation (13) being the Lorentz condition 

dA(x,t) | i a$(x,t) _ Q 

dx c dt 

In order to obtain the wave equation for the "vector potential" A, we write, 



c dt 

and equation (14) becomes, 



ldaM =4nJ( X ,t) (17) 



d 2 A(x,t) 4vr 

ax 2 c 2 «% 2 " c l ' j 1 j 



and similarly writing 

1 da(x, t) 



-Anp(x,t) (19) 



c 5x 

equation (15) becomes the wave equation for scalar potential $(x, t), 

d 2 $(x,t) I d 2 $(x,t) 



dx 2 c 2 dt 2 
The two definitions (17) and (19) imply that 



4irp(x,t). (20) 



dJ(x,t) dp(x,t) _ 

dx + dt " U {ZL) 

which is the equation of continuity. These considerations may be generalized to three space 
dimensions. 

The objective of the above long review is to stress the interesting fact that many of the 
fundamental linear equations of physics are obtainable from an elementary consideration of 
stochastic process. Of course, by no stretch of imagination, we may expect the whole of 
physics to follow from such a consideration. 



III. NONLINEAR EQUATIONS 



Nonlinear partial differential equations appear in all branches of physics and some of 
these have interesting properties such as soliton like solutions, infinite number of conserved 
objects and so on. The nonlinear diffusion equation in the form (in our notation) 

dP{x,t) 



dP(x,t) _ d 
dt dx 



f(P> 



dx 



(22) 



is well known in literature [|6|-[T0[| and the properties of its solutions have been extensively 
studied. 

Now we proceed with an aim at getting the above nonlinear equation and others out of the 
master equation (1), by suitable modifications. If we consider this to be a phenomenological 



equation, without any recourse to Poisson's distribution, then the obvious way to introduce 
nonlinearity is to introduce functions of x and t as multiplicative coefficients on right hand 
side of equation (1). As an example, suppose we apply randomness to Ax itself. We know 
that if a is chosen at random in the interval (0, 1) then the probability of a n to be in the 
interval (x, x + dx) is given by ^ 1+ ^ log2 f° r large n, and hence a possibility is to use ( 1+ f^ og2 
instead of Ax in equation (5), which is a special case of equation (1). Perhaps it would be 
simplest to replace Ax by P(x,t)Ax in (5), and the resulting equation is 

(23) 

where 

D= ( ^. (24) 
At v ' 

Or else, we may treat x and t in the same footing, that is set At — > P(x,t)At and Ax — > 
P(x, t) Ax instead of only Ax — > P(x, t)Ax and we get, 

^>, DPM ^H, (25) 

both equations (23) and (25) being nonlinear equations. Henceforth we set D — 1. 

It should be noted that this does not mean any nonlinear equation can be obtained this 
way. The condition that for At = , Ax = both left and right side of equation (1) must 
match is quite a restriction. However, by making use of equation (7) one may get many 
more equations by setting the source term as function of x, P and its derivatives or their 
combinations. This would be analogous to adding terms to the Lagrangian arbitrarily in the 
conventional method of getting equations of motion. 

It is an interesting fact that if the master equation (1) is modified in the first term of 
right hand side in the following way 

P±(x,t + At) = P±(x^Ax,t)(l - P + At) +P T (x± Ax,t)aAt (26) 

we get a nonlinear form of the Dirac equation in one space dimension, 

^ = . P} . v ^ +aP _ (27) 

dP- dP- , s 

— = - P+P _ + aP++v — (28) 



and by iteration a nonlinear analogue of telegraphers' equation results, 

dt* V dx* ~ F+ + F+ dt +VF+ dx +aF+ (29) 

Further generalizations would be to consider P as a complex multi- component object 
and readers may amuse themselves by putting objects such as supersymmetric variables, 
Pauli and other matrices etc. as coefficients of Ax in equation (1). 



IV. GROUP ANALYSIS 

Before we analyze the nonlinear equation (23) or (25) we would like to emphasize the 
important developments regarding the linear diffusion equation (7). It is well known that 
this equation goes over to Burgers equation under the Cole-Hopf transformation. Burgers 
equation and similar integrable ones have been much studied in recent decades for their 
importance in physical sciences for the existence of soliton like behaviour, infinite num- 
ber of conservation laws as well as for their deep mathematical connections to Backlund 
transformations, conformal invariance and so on. 

Equations of the form (23) have been analysed by Munier etal || and by Hill || in detail. 
It is found that the nonlinear diffusion equation of the form 



dip l2 d 2 i> 
~dt = ^dP 



^ (30) 



are equivalent to the classical diffusion equation for P, 

dP d 2 P 



dt dx 2 



(31) 



if we introduce x such that 



(32) 



and every nonlinear diffusion equation of the form 

(33) 

can be transformed to the following equation with a simpler nonlinearity 



dP _ d 
dt dx 



dP 



where ip(P,t) is the flux associated with equation (33). Hence, for this special case the 
analysis would be similar to that of the linear diffusion equation. 

However, in general the simplest nonlinear equation we would get from the master equa- 
tion by replacing Aa; — > f(P)Ax in equation (5) would be, 



as in equation (23) or 



BP 8 2 P 



dP d 2 P 



as in equation (25). 

Now we proceed to analyse the properties of the solutions of (35) by means of one param- 
eter groups as in [|BHlO|l. For the single dependent variable P and for the two independent 
variables x and t we have one parameter groups of the form, 

X! = f{x, t,P,e) = x + e£(x, t, P) + 0(e 2 ), 
h = g(x, t,P,e)=t + erj(x, t, P) + 0(e 2 ), 

Pi = h(x, t,P,e) = P + e((x, t, P) + 0(e 2 ). (37) 

We follow the standard procedure [P,|l(J to obtain the similarity variable and functional form 
of the solution by solving the first order partial differential equation 

dP dP 
t, P)-^ + V (x, t,P)-^= C(x, t, P) (38) 

for known functions £(x, t, P), T}(x, t, P) and ((x, t, P). Let 

d d d 

v = *• p) d^ + r]{x ' t ' p) di + c(x ' *' p) dP (39) 

be a vector field on the space X x , where coordinates represent the independent variables, 
the dependent variables and the derivatives of the dependent variables up to order 2. All 
possible coefficient functions £, rj, ( are to be determined so that the one parameter group 
exp(ev), thus obtained would be the symmetry group of the nonlinear equations (35) for the 
diffusion case and (29) for telegrapher's case. 

The determining equation for the symmetry group for the diffusion with non linearity, 
equation (35) is 



monomial 



coefficients 



a 2 p dp 

dxdt dx 


rjp = U 


( A ) 


d 2 P 
dxdt 


^ = u 


( B ) 


(dP\ 3 
\ dx I 


4pp - U 




(dP\ 2 dP 

\dx ) at 


t^pp = 


( D ) 


(dP\ 2 
\ dx > 


/a ot i n 

Isp — ^Sajp — U 






+ Vp — o 




(dP)(9P\ 

v dx A at ) 


-£p = -2Vxpf(P) ~ 3^p 


(G) 


ap 

dx 


& /(/M(2c,,.-C,! 


(H) 


dP 
dt 




(I) 


P° 


C* - = o 


(J) 



where prime denotes differentiation with respect to the argument and subscript denote 
differentiation with respect to the indicated variable. These equations turn out to be the 
same as those of the nonlinear diffusion equation of the form, 



dP _ d 
dt dx 



dP' 



considered in Hill |J. 

From (A), (B), and (G) it is easily seen that 



(40) 



£ = £0M), 



v = v(t) 



(41) 



Cp = 2£, + r, 



(42) 



where r is a constant. So 



Cpp = 0. 



(43) 



From (/) we get, 



(44) 



so that either 



or 



' f(py 
f(p). 



(46) 



pp 



that is, 



f(P) = a(P + b)' 



(47) 



where a, b and m denote arbitrary constants. If equation (45) holds, then from (h) and 
equation (44), we obtain, 



£(x,t,P) = P + >yx 
v(x,t,P) = 26 + 2-ft 
C(x,t,P) = 



(48) 



where, /3, 9 and 7 are arbitrary constants. 

Hence, the Lie algebra of infinitesimal symmetries of the equation is spanned by the 
three vector fields, 





d 


Vi = 






d 


v 2 = 


w 




d 


v 3 = 


dx 



f 2t 

and the commutation relations are given by, 

[vi,v 2 ] = 0, [Vi,V 3 ]=Vi, 



d_ 

di' 



(49) 



[v 2 , v 3 ] = 2v 2 . 



(50) 



The one parameter groups G,i generated by the Vj are given below. The entries give the 
transformed points exp(evj)(x, t, P) = (xi,ti,Pi). 

d : (x + e,t,P), 
G 2 : (x,t + e,P), 

G 3 : (e e x,e 2e t,P). (51) 



Each group Gi is a symmetry group and if P = q(x,t) is a solution of our nonlinear 
diffusion equation, so are the functions, 



pW = q( x - e,t), 
p( 2 ) = q( x ,t-e), 
p(3) = q (e- e x, e~ 2( 



t). 



(52) 



The groups we obtain are the same as those for equations (40) and so is the similarity 
variable M 
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x + a 



However, the functional form 



(53) 



p = sM 



of the solution satisfies the ordinary differential equation 

n , d 2 s ds 

whereas that corresponding to equation (40) is given by, 

aur as \ aa; / duj 



In the case /(P) is given by equation (46), 

m 

and for the time derivative of £ we get, 



2 <9£ _ cfy 



/(P) 



n? 



dx 



2 ' 



(54) 



(55) 



(56) 



while substituting (57) into (J) and using (58) gives, 

m 

So there are two possibilities arising out of equation (58), either for all constants m 

d£ d 2 i d 2 rj 



or for m = 4 , 



dt dx 2 dt 2 



d£ <9 3 £ d 2 r] 
dt dx 3 dt 2 



0, 



0. 



(57) 



(58) 



(59) 



(60) 



(61) 



Thus for all m we have, 

•q{x,t,P) — v + pt, 

C(x t t,P)=^^j{2a-p), (62) 

where p, v, a and p are arbitrary constants and the infinitesimal symmetries are spanned 
by four vector fields 



d_ 

dx 
d 

V2 = ^' 



- x 9 2 (P b) 9 
Va X dx m dP' 

d (p + b) d 

V4 = TPS'- 63 

at m oP 

and the commutation relations are given by, 

[vi, v 2 ] = [vi, v 4 ] = [v 2 , v 3 ] = [v 3 , v 4 ] = 0, [vi, v 3 ] = vi, [v 2 , v 4 ] = v 2 . (64) 

The one parameter groups Gi generated by the Vj are, 

G l : (x + e,t,P), 
G 2 : (x,t + e,P), 
G 3 : (e € x,t,(P + b)e%), 

G 4 : (x,e et ,(P + b)e'^), (65) 

and if P = y(x, t) is a solution to our non linear diffusion equation, so are the functions 

pW= y (x-e,t,P), 
pW=y(x,t-e,P), 
p( 3 ) =y( e - t x,t, (P-b)e~™), 

p( 4 ) = y( x ,e-%(P-b)e™). (66) 
The similarity variable in this case is given by, 



and the functional form of the solution is, 



P=U+-\ s{u) - b. (68) 

Now for the nonlinear form of the telegrapher's equation (29), arising out of the nonlinear 
Dirac equation (27), the independent determining equations of the symmetry group are given 
below. 



monomial coefficient 



d 2 P+ dP+ 
dt 2 dt 

d 2 P+ dP+ 

dxdt dt 

d 2 P+ 

dxdt 

d 2 P+ 
dt 2 


VP+ = o 
£p + = o 

& = v 2 r] x 
Vt = Zx 


(a) 
(b) 
(c) 
(d) 


9P+\ 
dt > 


2C*p + - Vtt + v 2 V xx + P+((p + - 2^)- 

P + Vf] x - P + T] x = 


(e) 


9P+ 
dx 


ftt + v 2 (2( xP+ + C+ 

p +v (Cp + - + p + ((p + - Q+ 






vP4( P+ -tx) = o 


(f) 




(tt-v 2 ( xx -Pl(( P+ -2Z x )+ 
P + a 2 (( P+ - + 3P + 2 C + P+C x + 






v( x p+ + a 2 ( = 


(g) 



The solutions are given by, 

£(x,t,P+) =Av 2 t + B, 
r}(x,t,P + ) = Ax + E, 

{(x,t,P + )=0, (69) 

where, A, B, and E are arbitrary constants and the infinitesimal symmetries are spanned 
by the three vector fields 



vi = — — , space translation, 
ox 

v 2 = — , time translation, 
ot 

d d 

v 3 = v 2 t— — h x—, hyperbolic "rotation" in x, t space, (70) 
with the commutation relations, 

[vi, v 2 ] = 0, [vi, v 3 ] = v 2 , [v 2 , v 3 ] = v 2 v 1 (71) 
The one parameter groups Gi generated by the are, 

d : (x + e,t,P + ), 
G 2 : (x,t + e,P + ), 

G 3 : (x + v 2 et,t + ex,P+). (72) 
This implies that if P + = z(x, t) is a solution to equation (29), so are the functions 

P+i = z(x - e, t), 
P+2 = z(x,t- e), 

P +3 = z(x — v 2 et, t — ex), (73) 

where e is any real number. 

In order to compare the above vector fields of equation (70) with those of the linear 
second order form of the telegrapher's equation (3), we have the corresponging independent 
determining equations of the symmetry group: 



monomial 


coefficient 








d 2 P+ 8P+ 
dx 2 dx 




= 






(P) 


d 2 P + 8P+ 
dx 2 dt 


Vp + ~- 


= 






(Q) 


d 2 P+ 
dx 2 


tx = 


Vt 






(R) 


d 2 P+ 
dxdt 


e* = 


V 2 T)x 






(S) 




Cp+P- 


f - ZVtP 


H + 4an P+ = 




(T) 


0P+ 
dx 


Ctt — 


V 2 ixx + 


2v 2 ( xP+ + 2a( t 




0(U) 


dP+ 
dt 


Vtt - 


v 2 rj xx - 


2( t p+ - 2ar]t = 





(V) 


(P + )° 


Ctt — 


V 2 Cxx + 


2a( t = 




(W) 



The solutions are given by, 



£(x,t,P+) =Kv 2 t + L, 
r](x, t, P + ) = Kx + M, 

C(x, t, P + ) = -KaxP + + NP+, (74) 

where K,L,M, and N are arbitrary constants. The infinitesimal symmetries are spanned by 
the four vector fields 



1 dx ' 



V2 dt 



V 3 = P+ 



0_ 

)x 
d 
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d 



dP+ 



Wi = vh Tx +x Wt- axP+ ^ (75) 

with commutation relations 

[vi, v 2 ] = [vi, v 3 ] = 0, [vi,v 4 ] = v 2 - av 3 , 

[v 2 , v 3 ] = 0, [v 2 , v 4 ] = v\ u [v 3 , v 4 ] = 0. (76) 
We have ignored the obvious infinite dimensional subalgebras in the above analysis. 



V. CONCLUSION 



The main objective of this paper was to extend the method of deducing some fundamental 
linear partial differential equations of physics from a consideration of stochastic arguments to 
the nonlinear case. We saw that this could be achieved in a very simple way by modifying the 
master equation to obtain "nonlinear diffusion" equation, a "nonlinear Dirac equation" in 1 
+ 1 dimensions and the corresponding "non linear telegrapher's equation" . As a preliminary 
step towards the analysis of the properties of the solutions, we have considered the group 
classification problem of the first and the last one by means of one parameter groups. The 
infinitesimal symmetry group of the nonlinear telegrapher's equation is spanned by a vector 
field corresponding to a "hyperbolic rotation" of x and t. For our type of diffusion equation, 
though the group structure is similar to that of standard nonlinear diffusion equation, the 



ordinary differential equations obtained is different and the results are similar when m = 4 
in our case, but m = — | in the standard case (m being the highest power of the dependent 
variable in coefficient to the term in the nonlinear diffusion equation). The physical 
applications of this equation has been widely studied in the context of gas dynamics and 
plasma physics etc. We expect the other two equations to have similar important applications 
in physics with rich mathematical structure and we leave it for future study. However, as a 
comparison of equations (70) and (75) shows one does see explicitly which symmetries gets 
broken when the equation is modified. 
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